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The order of the chiral transition in A^y = 2 QCD* 
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The order of the chiral transition for Nf = 2 is an interesting probe of the QCD vacuum. A strategy is developed 
to investigate the order of the transition using finite size scaling and its relation to color confinement. An in-depth 
numerical investigation has been performed with KS fermions on lattices with Nt = 4 and Ns — 12, 16, 20, 24, 32 
and quark masses am, ranging from 0.01335 to 0.35. The specific heat and a number of susceptibilities have 
been measured and compared with the expectation of an 0(4) second order and a first order phase transition. A 
second order 0(4) is excluded, whilst data are consistent with a first order. 



QCD with Nf = 2, is a key system to under- 
stand confinement. The phase diagram for two 
flavors of mass degenerate quarks is schematically 
depicted in Fig. ^ where the deconfinement tem- 
perature Tc is plotted versus the quark mass and 
the baryonic chemical potential fi. At /i = the 
phase transition is well understood at high masses 
(m > 2.5 GeV), where the quarks decouple: as in 
the quenched limit m ^ oo the transition is first 
order and (L), the Polyakov line, is a good order 
parameter. At rn ~ a chiral transition exists, 
where chiral symmetry is restored, and the chiral 
condensate (ipip) is an order parameter. At some 
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Figure 1. Schematic phase diagram of Nf — 2 
QCD. 



temperature also the axial C/a(1) is expected to 
be restored, and in fact the topological suscepti- 
bility drops to zero around Tc ^ . In principle at 
m ~ there are 3 transitions (chiral, axial [/^(l), 
deconfinement): it is not clear if they coincide, 
and the question cannot even be asked if there 
is no independent definition of deconfinenement. 
An effective description of the chiral transition 
can be given in terms of an effective free energy, 
by usual renormalization group plus e-expansion 
techniques |3- Under the assumption that the 
scalar and pseudoscalar modes are the relevant 
critical degrees of freedom, there exist no infrared 
stable fixed points for Nf > 3, so that the transi- 
tion is expected to be first order. For Nf = 2 the 
transition is first order if the anomaly is negligible 
(rrir,' w 0) at Tc (the symmetry is 0(4) x 0(2)); 
it can be second order with symmetry 0(4) if the 
anomaly survives the chiral transition. In the first 
case the transition surface around m = is first 
order, and there exists no tricritical point in the 
/i — T plane. If the chiral transition is instead sec- 
ond order, then by general arguments, the surface 
is a crossover, and a tricritical point is expected 
in the fJ. — T plane (see e.g. 0)- 

If the latter is the case, the decoirfining tran- 
sition cannot be order-disorder, there is no order 
parameter for confinement, and a state of a free 
quark can continuously be transfered below the 
"deconfining temperature" . The transition line 
in Fig. n is defined by the maxima of a num- 
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Figure 2. Specific heat (top) and Xm (bottom) peak value for Runl (left) and for Run2 (right), divided 
by the appropriate powers of Lg (Eqs. 113114(1 to give a constant. Horizontal stripes indicate the Icr 
confidence region for a fit with a constant value, x^/do/ is also shown. 



ber of susceptibilities (the specific heat Cy, the 
chiral susceptibilitiy, . . . ) which happen to co- 
incide within errors. Existing literature on the 
subject is admittedly not conclusive about the 
order of the chiral transition | 4l5l6l7l5| . with a 
diffuse tendency to assume a second order 0(4) 
and crossover. The problem is fundamental and 
deserves additional study. 

1. Strategies 

The order and the universality class of the tran- 
sition can be determined by a finite size scaling 
analysis of susceptibilities. The approch is based 
on the renormalization group, and is valid if the 
correlation length of the order parameter ^ goes 
large with respect to the lattice spacing a, so that 
a/^ « 0: this is true for second order and for weak 
first order transition j9lim . A key quantity is the 
specific heat Cv, for which 

Cy - Co = L^/''$c(Tiy^ am,Lf ) (1) 



T = 1 — T /Tc is the reduced temperature, the 
index v is defined by the critical behavior of ^, 
^ oc |t|~'', and Co is a subtraction constant due to 
additive renormalization ^H!- Eq. ^is valid inde- 
pendent of the knowledge of the order parameter. 
It must be stressed that the order of the transi- 
tion can in priciple be established regardless of 
any prejudice by looking at the behavior of the 
specific heat, which is therefore be the quantity 
of reference. 

If the order parameter is known the following 
scaling law is expected for its susceptibility 

X^Ll'^^^{TLl/-,am,Ly-) (2) 

V, a, 7 identify the order and universality class 
of the transition. A good order parameter must 
obey Eq. |21 

On the lattice 1/T — Lta{P,m), so that 

^=l-#4«^[/3e-/3 + M (3) 
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Figure 3. Peak of the specific heat divided by the appropriate power of the mass (Eq. IT^ to give a 
constant, for different scaling hypotheses. 



for sufficiently small masses while higher order 
terms in m can enter at higher masses. Here 



C = 



dhia{l3, m) 



d(3 





d\na 


, k = 


dm 


5 In a 


■J(/9,0) 


L a/3 J 



(4) 



In the quenched case k — and t cx /3c — /3, as 
usually assumed. In the existing literature for 
Nf = 2 the m dependence in Eq. is usually 
neglected, and r is set proportional to /3c — (3. 

Our strategy has been to assume the 0(4), 
0(2) exponent yh = 2.49 (see Table^ and to run 
by choosing arriq, such that anigL^'^ = const. 
The problem has then only one relevant scale and 
v, a, 7 can be determined: if 0(4) or 0(2) is the 
correct universality class the peaks should scale 

as (Cv - Co)f™'= oc L"^", x^'"'' <^ ■ More- 
over, by analiticity arguments |S| in the m — > 
limit one expects, if the transition is second order, 
that 

Cv-Cn = (amq)-'''^''y-^<^'c{TLl''') (5) 
Xrn-Xo = {am,r-^'^''y-^-^'^(TL\''') (6) 

The pseudocritical coupling should then obey the 
scaling law 



/3c - /3pc + km~ k'/Ll'" 



(7) 



In the case of a first order transition the 
analiticity arguments used in deriving the scaling 
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Critical exponents. 





Vt 


Vh 
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0(4) 


1.34 


2.49 


0.75 


-0.23 


1.48 


0(2) 


1.49 


2.49 


0.67 


-0.01 


1.33 


MF 


3/2 


9/4 


2/3 
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V^Order 
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3 


1/3 
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Eas. 15161 are no longer valid, however they are ex- 
pected to hold for small volumes compared to the 
critical volume. For larger volumes the suscepti- 
bility are expected to increase with the volumes. 

Section contains the results of the analysis. 

An additional susceptibility has been ana- 
lyzed PJj, i.e. p = 91n(/i)/9/3, where (/x) is an 
order parameter for confinement detecting dual 
superconductivity of the vacuum. 

2. Monte Carlo Simulations 

The lattice action used in Monte Carlo simu- 
lations was the standard staggered action. Con- 
figurations were updated using the Hybrid R al- 
gorithm and taking care of systematic errors due 
to finite molecular dynamics integration step size 
and finite conjugate gradient inversion residue. 
Two different sets of MC simulations were per- 
formed, called in the following Runl and Run2, 
having fixed for each set the value of aniqLy'^ as 
explained above while we set Lt = 4. The pa- 
rameters of our runs and the correspondig total 
number of trajectories collected are reported in 
Table 121 For each value of am^ and a number 
(10-^20) of simulations at different /3 values were 
performed in order to inspect and to have un- 
der control the whole interesting critical region. 
The multi-histogram reweigthing technique was 
exploited to obtain information at intermediate /3 
values. The following quantities were measured: 



Cv = 

Xm = 

Xij = 

Xee = 



1 

T 52 



InZ 



Xij I Xee ; Xie 



(8) 



V drriq 



lnZ^V[{{y^ijy)-{^^y] (9) 



VmP,) - {P,) {P,)], i,j^a,r (10) 
Vid^Do^ f) - i^Do^)^] (11) 



4 




X^e = V[{n{i^Doi^)) - {P^ {i^Dai^)] (12) 

where V = L^Nt is the volume; Dq is temporal 
component of the Dirac operator; Po-, Pt indi- 
cate the average spatial and temporal plaquette 
respectively; Cy denotes the specific heat which 
is a function of Xy, Xee and Xie! Xm is the sus- 
ceptibility of the chiral condensate. Only the dis- 
connected component of Xm is considered in the 
present work since it gives the dominant contri- 
bution as it results from previous works and as it 
was also checked for a fraction of our MC simu- 
lations. 

3. Numerical Results 

Having fixed anigLy'^ — const with yh = 2.49 
that is the value expected for 0(4) and 0{2) 
critical behavior, the following scaling formulas 
should hold (see Eas. lll2|l : 

Cv{T,L)-Co = L'^/-'<^c{rLl/n (13) 
Xm{r,L)-xo^L-y/''^^{TLl^n (14) 



In particular for the peaks of the susceptibilities 
we expect: C^^"'' oc L"/'' and x^"'' « Lt/". This 
behavior is only expected for the singular part of 
these quantities as the subtraction of constants 
Co and xo indicate in Eas. 1131 HI The measured 
peak values from data reweightening for the spe- 
cific heat and Xm for Runl and Run2 are shown 
in Fig. 13 The 0(4) and 0(2) critical behavior is 
clearly in contradiction with the observed quan- 
tities. In particular 0(4) and 0(2) scaling pre- 



Table 2 

Run parameters for the numerical simulations. 



Ls ■ TOjT varies in the range 


8.9,15.8] 




aniq 


# Traj. 


Ls 


Runl 


Run2 


Runl 


Run2 


12 


0.153518 


0.307036 


22500 


25000 


16 


0.075 


0.15 


87700 


131390 


20 


0.04303 


0.08606 


14520 


16100 


32 


0.01335 


0.0267 


14500 


15100 
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Figure 5. Determination of the reduced temperature for 0(4) and 1^* order. The smaU m region is 
enlarged in the inset picture. 



diets no singular behavior in the L —t oo limit for 
the specific heat as a is negative. Also for Xm 
the predicted exponents fail to reproduce lattice 
data. In the whole study a non-improved action 
was used. One could then wonder whether this 
disagreement is a consequence of a poor conver- 
gence to the continuum limit. However it seems 
unlikely that an infrared qualitative feature of the 
theory, like the growth of the specific heat, can 
heavily rely on improvement of the action. 

Wishing to test scaling behaviors other than 
the two previously considered, namely first or- 
der and mean field, we considered the scaling 
relations Eqs. 15161 valid in the limit aniq 0. 
Smaller lattices and larger masses were discarded 
for this analysis in order to satisfy the constrains 
of above equations and an additional set of sim- 
ulations with Ls = 24 and aruq = 0.04444 was 
added. We repeated the analysis of the scaling 
of the peaks of susceptibilities and the results are 
shown in Fig. El for the specific heat. Xm shows 
a similar behavior. From the value of the x^ /dof 
is easly seen that the only choice compatible with 
the data is a order scaling. Although the up- 
per mass limit of validity of Eas. l5l6l is not known 
a priori, it should be noted that restricting the 
mass range of the fit worsen the agreement with 
0(4), 0{2) and MF hypotheses while for first or- 
der it remains unchanged. The scaling of suscep- 
tibilities at all fi can also be investigated. Fig. ^ 
shows the quality of the scaling for the two dif- 
ferent working assumptions of 0(4) and I''* order 
both for the specific heat and Xm- Even if it is 



difficult to make precise quantitative statements 
about the quality of these scalings, the first or- 
der case seems to give better results for both of 
the susceptibilities considered. This is clearly the 
case for the height of the susceptibilities as shown 
above, but it is also true for the width of these 
curves. The horizontal scale of these plots is given 
by the scaling variable tL}J^ where r, the reduced 
temperature, depends on the unkwon value of k 
(see Eq. O. The actual value of k was extracted 
form the data with the fits shown in Fig.O Both 
0(4) and first order critical exponent allow a good 
fit of r so that it is not possible to discriminate 
between the two cases considering only pseudo- 
critical couplings. 

Having detected a signal of a first order phase 
transition, it is natural to look at time histo- 
ries of observables to find signals of metasta- 
bilities. Such metastabilities should manistifest 




Figure 6. MC histories and observable histograms 
showing a double peak structure. Plaquette is ac- 
tually shown and chiral condensate have a similar 
behavior. 
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themselves as double peaks structures in the his- 
tograms of observable distributions. Indeed we 
found such structures for some of our time his- 
tories at largest volumes and smallest masses, as 
shown in Fig. |B| Even though this is a prelim- 
inary indication of a first order phase transition 
it is not conclusive and we were not able to evi- 
dence this double peak structure consistently, in 
particular in the data taken at avriq — 0.0267 with 
Ls — 32, where it was expected to be present. 

4. Conclusions 

We have presented a systematic study address- 
ing the fundamental question of the order of 
the chiral transition in two-flavor QCD. A large 
numerical effort has been done with the aim 
of clarifing the nature of the chiral transition. 
Our MC simulations were run on APEmille ma- 
chines with a total computation cost of about 
5 TFlops X month. 

The hypothesis of a second order phase transi- 
tion belonging to the universality class of 0(4) or 
0(2) has been carefully tested on the lattice by 
looking at the finite size scaling of the basic ther- 
modynamic susceptibilities, first of all the specific 
heat and then the chiral condensate susceptibil- 
ity. To do that, we have fixed the scaling variable 
am^L^'' in Eas. lll2l and we have studied the resid- 
ual Lg dependence. Our results, presented in the 
first part of Sect. 13 seem to exclude both the 0(4) 
and 0(2) universality classes. In particular, the 
growth of the specific heat going towards small 
masses and large volumes is in constrast with the 
expected behavior of a 0(4) or 0(2) transition in 
which there should be no singular contribution to 
Cv- The finite size scaling of the chiral conden- 
sate susceptibility is also incompatible with 0(4) 
and 0(2) universality class. 

An extented analysis at small masses, pre- 
sented at the end of Sect. |31 clearly shows that 
a first order phase transition is instead consistent 
with lattice data while all other open choises are 
disfavored. Together with the scaling of the spe- 
cific heat and order parameter susceptibility, also 
the time hystories of observables show some sig- 
nals of metastabilities, i.e. double peak structures 
in distributions, which, although not conclusive, 



are a qualitative indication of a first order phase 
transition. 

We have used the standard staggered action 
and Lt = 4. We plan to repeat the analy- 
sis with larger Lt and improved action. How- 
ever the growth of the peaks which excludes 0(4) 
and 0(2) is a typical infrared phenomenon, most 
likely independent on short range improvements. 
We also plan to run keeping the scaling variable 
OTOgL^'' constant under the assumption of a first 
order transition, to check directly if it is consis- 
tent. The indication JI] coming from the dual 
superconductivity disorder parameter is also in 
the direction of a first order phase transition. 
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